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Abstract 

The determination of the magnitude and "sign" of the J/ip polarization crucially 
depends on the reference frame used in the analysis of the data and a full under- 
standing of the polarization phenomenon requires measurements reported in two 
"orthogonal" frames, such as the Collins-Soper and helicity frames. Moreover, 
the azimuthal anisotropy can be, in certain frames, as significant as the polar 
one. The seemingly contradictory J/ip polarization results reported by E866, 
HERA-B and CDF can be consistently described assuming that the most suit- 
able axis for the measurement is along the direction of the relative motion of 
the colliding partons, and that directly produced J/V>'s are longitudinally po- 
larized at low momentum and transversely polarized at high momentum. We 
make specific predictions that can be tested on existing CDF data and by LHC 
measurements, which should show a full transverse polarization for direct J ftp 
mesons of pt > 25 GeV/c. 



Submitted to Phys. Rev. Lett. 



^ Laboratorio de Instrumentagao e Fisica Experimental de Partfculas (LIP), 
Lisbon, Portugal 

2 ) CERN, Geneva, Switzerland 

3 ) Instituto Superior Tecnico (1ST) and Centro de Fisica Teorica de Partfculas 
(CFTP), Lisbon, Portugal 



The existing measurements of J/tp polarization in hadronic collisions represent 
one of the most difficult challenges currently faced by models of quarkonium produc- 
tion (see, for example, Refs. [U [2] and references therein). The often emphasized 
disagreement between experiment and theory is, however, only one aspect of the 
problem. The experimental knowledge itself looks contradictory when different po- 
larization measurements are compared, in terms of "sign" , magnitude and kinematic 
dependence, as illustrated in Fig. CQ, which shows the data reported by CDF [2], 
HERA-B [3] and E866 [TJ. 
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Figure 1: A# versus pr, as reported by E866, HERA-B and CDF (statistical and 
systematic errors added in quadrature). 

Besides the obvious interest of understanding the mechanisms of quarkonium po- 
larization, having a clear (data driven) description of the polarization measurements 
is also important to evaluate detector specific corrections needed to extract physics 
results from the data. Production cross sections, for instance, might significantly 
depend on the polarization scenario used in the calculation of acceptance corrections. 
The polarization measurements are undeniably complex and involve difficult exper- 
imental problems. There is, however, an additional cause for the blurred picture 
emerging from the comparison of the existing measurements: different experiments 
have often chosen different polarization frames to perform their analyses. The in- 
fluence of such choices on the measured angular distribution of the decay leptons 
is generally underestimated. In fact, different analyses of the same two-body an- 
gular decay distribution may give qualitatively and quantitatively different results 
depending on the definition of the polarization frame. 

Several polarization frame definitions have been used in the past. In the helicity 
frame the polar axis coincides with the flight direction of the J/ip in the centre-of-mass 
frame of the colliding hadrons. A very different approach is implicit in the definition 
of the Collins-Soper [5] frame, where the polar axis reflects, on average, the direction 
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of the relative velocity of the colliding partons, the approximation being especially 
good if we can neglect the smearing effect due to the parton intrinsic transverse 
momentum. We denote by d the angle between the direction of the positive lepton 
and the chosen polar axis, and by (p the azimuthal angle, measured with respect to 
the plane formed by the momenta of the colliding hadrons in the J/i/j rest frame (the 
"production plane"). The angular decay distribution, symmetric with respect to the 
production plane and invariant under parity transformation [3, [6] , is usually defined 
as: 



If the J/tp is observed in a given kinematic configuration, any two polarization frames 
differ only by a rotation around the axis perpendicular to the production plane (the 
u y axis"). The functional dependence of the decay distribution on the angles $ 
and <p is invariant with respect to such a rotation, but the numerical values of A#, 
A,^ and X v change in a correlated way. In particular, a rotation by the angle 5 = 
1/2 arctan[2 X^ v ,/(X cp — A#)] (or 45° when A^ = A#) leads to a frame where A^ is 
zero, i.e., a frame with axes along the principal symmetry axes of the polarized 
angular distribution. The experimental determination of \# v can, therefore, provide a 
criterion for the choice of a particularly convenient reference frame for the description 
of the angular distribution. 

While all three coefficients provide interesting and independent information, most 
available measurements of J/-0 polarization are restricted to A#. This limits the 
possible interpretations of the results and forces us to rely on model-dependent as- 
sumptions when comparing results obtained by experiments using different reference 
frames. Even the seemingly simple classification of "transverse" or "longitudinal" 



dent when the decaying particle is produced with small longitudinal momentum, the 
Collins-Soper (CS) and helicity (H) polar axes becoming perpendicular to each other. 
In this case (assuming A^ = = 0, for simplicity), if in one frame we measure a 
value A,?, the value measured in the second frame is smaller and of opposite sign, 
A^ = — A#/(2 + A#), while an azimuthal anisotropy appears, A' = A#/(2 + A$). 

There is a further reason for performing the experimental analyses in more than 
one reference frame. The J/i/j acquires its polarization with respect to a "natural" 
polarization axis which is, a priori, unknown and not necessarily definable event- 
by-event in terms of observable quantities. In practice, a fine-grained scan of the 
multidimensional phase-space of the J/ip production process is not possible, due to the 
limited sample of collected events, which forces the decay distribution to be measured 
as an average over a wide spectrum of kinematic configurations. This means that the 
orientation of the polar axis of the chosen frame with respect to the "natural axis" 

1 Following a common (albeit misleading) practice, the polarization is defined as transverse (lon- 
gitudinal) when > (A^ < 0). 
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is, in fact, dependent on the reference frame. This is particularly evi- 
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changes from event to event, depending on the momentum of the produced J /if). The 
resulting superposition of many distributions, equal in shape but randomly rotated 
with respect to one another, is "smeared" into a more spherically symmetric shape. 
As a consequence, the measured absolute values of A# and A^ are smaller than what 
would be measured in a fixed kinematic configuration and in the "natural frame". 
Therefore, independently of any prior theoretical expectation, the frame closest to 
the natural frame is the one providing the smallest 5 angle and the most significant 
|A,?|. 

The HERA-B experiment recently reported the three coefficients determining the 
J /if) decay angular distribution, in three reference frames [3], providing a clear picture 
of how the shape of the distribution changes from frame to frame. Before discussing 
kinematical dependences, we start by considering the values integrated in the phase 
space window covered by HERA-B: in the CS frame, A# = —0.31 ± 0.05 and A^, = 
-0.02±0.02; in the H frame, A* = -0.11±0.05 and X v = -0.07±0.02 (statistical and 
systematic errors added in quadrature). Furthermore, 5 has a much larger error in the 
H frame (10° ±20°) than in the CS frame (3° ±3°), reflecting the poorer precision with 
which the "tilt" of a more spherically symmetric shape can be determined. With the 
largest \X.g\ and a A^ compatible with zero, the CS frame is shown by the HERA-B 
measurements to provide a simpler angular distribution than the H frame. 
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Figure 2: A.# as a function of p, from E866, HERA-B and CDF data (statistical and 
systematic errors added in quadrature). 

We now address the kinematical dependence of the J /if) polarization. Figured] 
shows that, in the CS frame, E866 pE] observed a small J /if) transverse polarization 
(A# ~ 0.1) while the HERA-B pattern indicates longitudinal polarization, of decreas- 
ing magnitude with increasing pr- These are not conflicting observations, given the 
significantly different Xp windows covered: the average J /if) longitudinal momentum, 
in the centre of mass of the collision system, is 7 and —1.4 GeV/c for E866 and 
HERA-B, respectively. Indeed, Fig. [5] shows that the total J /if) momentum (here cal- 
culated using average xp values) provides a good scaling between the two fixed-target 
data sets. As also shown in Fig. [TJ CDF [2] reported that, above p-p = 5 GeV/c, the 
J /if) polarization is longitudinal in the H frame, with A# steadily decreasing with 
Pt. To see how the CDF pattern compares to the fixed-target data sets, we need 
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to convert the published values to the CS frame. We did this translation (using the 
relations presented above) assuming that X v = in the CS frame, as suggested by the 
HERA-B measurements. The resulting pattern, seen in Fig. [2J is perfectly aligned 
with the HERA-B and E866 data points. 

This smooth overlap of the three data sets suggests a simple polarization scenario, 
where the CS frame is taken to be a good approximation of the natural polarization 
frame (A^ = 0, A^ = 0) and A# is a monotonically increasing function of the total 
J /ip momentum. Before searching a suitable function, we remind that a significant 
fraction of the observed J/ip mesons results from Xc and ip' feed-down decays [7]: 
/fd = 0.33 ± 0.05. Irrespectively of the possible polarizations of these charmonium 
states, it is reasonable to assume that the strong kinematical smearing induced by the 
varying decay kinematics reduces the observable polarization of the indirectly pro- 
duced J/tp mesons to a negligible level. The feed-down contribution from 6-hadron 
decays can also be neglected: very small at fixed-target energies and experimentally 
rejected in the CDF analysis. Therefore, the observed polarization should be essen- 
tially determined by the directly produced J/^'s. The curve in Fig. [2J represents a fit 
of all the data points using the simple parametrisation 
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where the polarization of the directly produced J/^'s changes from fully longitudinal 
at zero momentum to fully transverse at asymptotically high momentum. The fit 
gives p = 5.0 ± 0.4 GeV/c and k = 0.6 ± 0.1, with xV ndf = 3.6/13. 
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Figure 3: pt dependence of A# in the CS 
frame, according to Eq. [2jand as reported 
by E866 and HERA-B. 



Figure 4: p T dependence of A# in the H 
frame, as derived from Eq. [2J and as re- 
ported by HERA-B and CDF. 



Our simple parametrisation provides a good description of the existing data sets, 
as can be seen in Figs. [3] and HI where the widths of the bands correspond to ±1<t 
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variations in the two fitted parameters as well as in the J/ip feed-down fraction. The 
derivation of the A# pattern in the H frame (needed, in particular, to address the 
CDF case) incorporates the "kinematical smearing" induced by the decays and the 
differential acceptances of the experiments (using a simple Monte Carlo procedure). 
In the narrow rapidity window of CDF, where the maximum J ftp longitudinal mo- 
mentum (~4 GeV/c) is always smaller than the minimum j>t (5 GeV/c), the helicity 
and Collins-Soper frames are essentially orthogonal to each other. Therefore, the 
decrease of \# with pr seen in the H frame (Fig. H|) is equivalent to an increase in the 
CS frame, as shown in Fig. [5j 
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Figure 5: dependence of A# in the 
CS frame, calculated for the energy and 
rapidity windows of the PHENIX, CDF 
and CMS experiments. 
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Figure 6: Same as previous figure, but 
for A M in the H frame. 



Assuming that the decay distribution has a purely polar anisotropy in the Collins- 
Soper frame, with \$ depending on momentum according to Eq. [2j CDF should 
observe a significant azimuthal anisotropy in the helicity frame, with a A^ pattern 
(shown in Fig. [6]) similar in magnitude but of opposite sign with respect to their 
A,?(pt) curve. By simply repeating the J/ip polarization analysis using the CS frame 
and by reporting the azimuthal angular distribution, CDF can clarify whether the 
polarization of the J/xjj is, also at collider energies, induced along a direction close to 
the parton-parton interaction line. 

Figures [5] and [6] also show the calculated pr dependence of A#, in the CS frame, and 
of Ap, in the H frame, for the kinematical conditions of the PHENIX (y/s = 200 GeV, 
\tj\ < 0.35) and CMS (y/s = 14 TeV, \tj\ < 2.4) experiments. If Eq. [2] remains valid 
up to LHC energies, we should see A# saturating for p T values higher than those 
probed by CDF, with a magnitude determined by the fraction of directly produced 
J/-0 mesons. 
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We will now summarise our main messages. 1) To investigate the J/i/j polarization 
and understand its origin, it is essential to know both the polar and azimuthal distri- 
butions, and their kinematical dependences, in at least two frames. The Collins-Soper 
and helicity frames, exactly orthogonal to each other at mid-rapidity, represent a good 
minimal set of polarization frames. 2) The HERA-B measurements show a pure po- 
lar anisotropy in the CS frame while a mixture of polar and azimuthal anisotropies 
is seen in the H frame, indicating that the J/ip decay angular distribution assumes 
its simplest shape when observed with respect to a polar axis that reflects the rel- 
ative momentum of the colliding partons rather than the J/ip momentum. 3) The 
seemingly contradictory patterns published by E866, HERA-B and CDF can be con- 
sistently reproduced assuming that the polarization (in the CS frame) of the directly 
produced J /if) mesons changes gradually from fully longitudinal at zero momentum 
to fully transverse at very high p. 4) This suggests that the longitudinal polarization 
reported by CDF in the H frame is, in fact, the reflection of a transverse polar- 
ization (around twice as large) in the CS frame, increasing with p^. Moreover, an 
azimuthal anisotropy of the decay distribution should exist in the H frame, with the 
same significance as the polar result. 5) Our polarization scenario predicts that the 
polar anisotropy of the prompt J/ip sample will saturate, for px above ~25 GeV/c, 
at A,? ~ 0.6-0.7, a value determined by the magnitude of the if)' and Xc feed-down 
contributions, assumed to be of negligible observable polarization. This prediction, 
easily verifiable at the LHC, can be placed on more robust grounds once CDF re- 
ports the complete angular distribution in the CS frame or, at least, the azimuthal 
component in the H frame. 

We acknowledge very stimulating discussions with R. Spighi. P.F. and H.W. are 
supported by FCT (Portugal) contracts SFRH/BPD 42343/2007 and 42138/2007. 
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